a; 



Photometric and spectroscopic observations of the 
F3+M3 eclipsing binary T-LyrO-08070 

Omiir Qakirli'^'*, Cafer Ibanoglu*^, E. Sipahi*^ 

"•Ege University, Science Faculty, Astronomy and Space Sciences Dept., 35100 Bornova, 

Izmir, Turkey. 



Abstract 

The multi-color photometric and spectroscopic observations of the newly 
discovered eclipsing binary T-LyrO-08070 were obtained. The resultant light 
and radial velocities were analysed and the absolute parameters of the com- 
ponents were determined. The system is composed of an F3 and an M3 
main-sequence stars. Masses and radii were estimated to be 1.37±0.23 Mq 
and 1.60±0.09 Rq for the primary and 0.32±0.04 Mq and 0.86±0.06Rq for 
the secondary star. The less massive secondary component has a radius at 
least two times larger with respect to its mass. Using the BVJHK magni- 
tudes of the system we estimated an interstellar reddening of 0.22 mag and 
a distance to the system as 479±36 pc. 
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1. Introduction 



Transit events for the T-LyrO-08070 (Vmax=12.15 mag.) were detected 
during r egular operation s of th e Trans- Atlantic Exoplanet Survey network 



{TrES). 



Fernandez et al. 



(120091 ) classified this star long-period single- 
lined eclipsing binary with an F-star primary and an unseen M-dwarf sec- 
ondary. They obtained follow-up spectroscopic observations using 1.5 m 
Wyeth Reflector at the Oak Ridge Observatory and 1.5 m Tilinghast Re- 
flector at Whipple Observatory. In order to reveal fundamental character- 
istics of the T-LyrO-08070 they combined photometric data from TrES sur- 
vey and KeplerCam photometry with those single-order echelle spectra in a 
wavelength window of 45 A centered at 5187 A. They obtained the radial 
velocities only for the F-star. Since their spectra cover a narrow wavelength 
interval, there is a clear dependence between the effective temperature and 
the metal abundance of the star. For this reason they have, modeled the 
spectra for metallicity index of -1.0, -0.5, 0.0 and +0.5 to obtain effective 
temp erature, projected rotat ional velocity and surface gravity of the primary 



star. 



Fernandez et al. 



(I2OO9I ) r nodeled the Kep l erCam primary transit light 
curve employing the method of iMandel fc Agoll (120021 ). which was an analyt- 
ical light curve analysis developed for planetary transit searches. Assuming 
orbit-rotation synchronization and adopting Fe/H=-0.5, the mass and radius 
of the primary star were derived as 0.95 MQand 1.36 Rq. Under the same 
assumptions, they derived the mass and radius of the secondary star as 0.345 
Moand 0.360 Rq. 

We planned new spectroscopic and photometric observations of the T- 
LyrO-08070 to obtain the masses, radii and effective temperatures of the 
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primary and secondary stars. The photometric hght curves were obtained 
in the BVRI passbands and analysis are presented. The radial velocities for 
both components are also obtained and analysed. 



2. Data acquisition 

2.1. Spectroscopy 

Optical spectroscopic observations of the T-LyrO-08070 were obtained 
with the Turkish Faint Object Spectrograph Camera (TFOSC)|^ attached 
to the 1.5 m telescope on 8 nights (July, 2010) under good seeing condi- 
tions. Further details on the telescope and the spectrograph can be found at 



http://www.tug.tubitak.gov.tr;. The wavelength coverage of each spectrum 
was 4000-9000 A in 12 orders, with a resolving power of A/AA 7000 at 6563 
A and an average signal-to-noise ratio (S/N) was ~120. The exposure times 
were between 10 and 60 minutes during observations. We also obtained a 
high S/N spectrum of the M dwarf GJ 740 (MO V) an d GJ623 (Ml. 5 V) fo r 



use as templates in derivation of the radial velocities (iNidever et al 



20021 ) 



The electronic bias was removed from each image and we used the 'crre- 
ject' option for cosmic ray removal. Thus, the resulting spectra were largely 
cleaned from the cosmic rays. The echelle spectra were extracted and wave- 



length ca' 



ibrated by using Fe-A r lamp source with help of the IRAF echelle 



package (ITonry fc David . Il979l ) . 



The stability of the instrument was checked by cross correlating the spec- 
tra of the standard star against each other using the FXCOR task in IRAF. 



^http:/ /tug. tug. tubitak.gov.tr/rttl50_tfosc.php 
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The standard deviation of the differences between the velocities measured us- 
ing FXCOR and the velocities in 



Nidever et al. 



fl2002[ ) was about 1.1 km s~\ 



2.1.1. Spectral classification 

We have used our spectra to reveal the spectral type of the primary 
component of T-LyrO-08070. For this purpose we have degraded the spectral 
resolution from 7000 to 3 000, by convolving them with a Gaussian kernel of 
the appropriate width, and we have measured the equivalent widths {EW) of 
photospheric abs orption lines for the spect ral classification. We have followed 



the procedures of 



Hernandez et al. 



(120041 ) , choosing helium lines in the blue- 
wavelength region, where the contribution of the secondary component to the 
observed spectrum is almost negligible. From several spectra we measured 

1 A. From the 



WHei+FeiA4922 = 0.38 ± 0.02 A and EW,^^^s^^_^_b^b±_^ 



calibration relations i^VT-Spectral-type of iHernandez et al. I (120041 ). we have 
derived a spectral type of F3±2 for the primary star. In Fig.l we compare 
the spectrum of T-LyrO-08070 with the spectra of standard stars HD 224639 
(FOV), HD 223623 (F2V) and HD 18215 (F5V). The spectra of standard stars 
are taken from the ELODIE archive! database. This comparison confirms the 
spectral-type of the more massive star derived from the EW measurements. 

2.1.2. Rotational velocity 

The width of the cross-correlation profile is a good tool for the measure- 
ment of V sin i (see, e.g., Queloz et al. 1998). The rotational velocities {v sin i) 
of the two components were obtained by measuring the FWHM of the CCF 



•^http://atlas. obs-hp.fr/elodie/ 
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Figure 1: Comparison of the observed spectrum of the binary with those of the F(0-5)V 
spectral type standards in the AA4800 — 5300 A region. 



5 




5140 fjlf)!) 5160 5170 5160 

Wavelength (A) 



190 5200 



Figure 2: Observed spectrum (points) of T-LyrO-08070 in the vicinity of Mgl triplet. The 
synthetic spectrum (F3 V+M3 V) is displayed with continuous line. 

peaks in nine high-S/N spectra of T-LyrO-08070 acquired close to the quadra- 
tures, where the spectral lines have the largest Doppler-shifts. In order to 
construct a calibration curve FWHM-f sini, we have used an average spec- 
trum of the HD 128429, acquired with the same instrumentation. Since the 



rotational ve 
km s~^, e.g.. 



ocity of HP 



Rover et al 



2842 9 is very low but not zero (F5V, vsini ~14 



( 120021 )). it could be considered as a useful template 
for F-type stars rotating faster than vsini ^ 10 km s^^. The spectrum of 
HD 128429 was synthetically broadened by convolution with rotational pro- 
files of increasing v sin i in steps of 5 km s~^ and the cross-correlation with 
the original one was performed at each step. The FWHM of the CCF peak 
was measured and the FWHM-f sin i calibration was established. The v sin i 
values of the components of T-LyrO-08070 were derived from the FWHM 
of their CCF peaks and the aforementioned calibration relations, for a few 
wavelength regions for the best spectra. This gave values of 58±1 km s~^ 
for the primary and 25±9 km s~^ for the secondary star. The rotational 
velocity of the primary star is in a good agreement with that estimated by 



Hernandez et al. 



20041 ). 
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For construction of the synthetic spectrum of the system the spectra of 
the standard stars {l Psc and GJ740), obtained with the same instrumenta- 
tion, have been rotationally broadened by convolution with the appropriate 
rotational profile and then co-added, properly weighted by using physical 
parameters (Ti, T2, -Ri, R2, vsinii2) of the components as input parameters 
and, then, Doppler-shifted according to the radial velocities of the compo- 
nents. Ti and vsinii^2 were determined in §2.1.1 and §2.1.2. T2 and Ri^2 
were obtained from the combination of the light and radial velocity curves 
analyses (see §4). In Fig. 2 we compare the observed spectrum of T-LyrO- 
08070 around the Mg I triplet with the synthetic one for a system consisting 
of F3V+M3V type stars. 

2.2. Photometric observations 

In order to provide high-quality light curves of the system in different 
passbands, we used SI 1100 CCD Cameralfl mounted on the 1 m telescope at the 
TUBITAK National Observatory, Turkey. SI 1 100 Camera utilizes a Cryo- 
genic cooler 4kx4k, Fairchild 486 CCD that gives a 21' x 21' field and a 
pixel size of 0".72 when the binning is 2x2. Differential aperture photome- 
try was performed to obtain the light curves. The observed magnitudes of 
T-LyrO-08070 in the B and V passbands were compared with those for the 
comparison and check stars. The coordinates and magnitudes of the stars 
observed are given in Table 1 . The B and V magnitudes of the variable are 
obtained in this study. The typical standard deviations of the differential 
magnitudes are about 0.005 mag. 



^http://www. specinst.com/ 
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Figure 3: The residuals between the observed times of minimum hght and computed with 
the new ephemeris. The square refers to the secondary minimum. 

Table 1: The coordinates and BVJHK magnitudes of the stars observed. 
Star a (5 B (mag) V (mag) J (mag) H (mag) K (mag) 

T-LyrO-08070 19'' 19" 04" 38° 40' 57" 12.747 12.150 11.507 11.321 11.261 
GSC 3121-01011 19'' 19'" 50" 38° 40' 57" 10.81 9.92 8.173 7.694 7.633 
GSC 3121-00379 19'' 19'" 33 " 38° 40' 57" 11.7 11.5 10.703 10.569 10.544 



0.01 

S" 0.00 r 
I - 

Q 0.00 — 



-0.00 - 



2. 3. Effective temperature of the primary star 

We find an observed color of B-V=0.60±0.02 mag for the binary system 
T-LyrO-08070 at out-of-eclipse which should be reddened if the spectral type 
F3 is considered. The observed infrared colors of J-H=0.186±0.039 and H- 



K=0.060±0.035 given in the 2MASS catalog f ICutri et al. I . l2003[ ) correspond 
to a spectral type of F3±2 is in a good agreement with that we derived by 
spectroscopy. 

The effective temperature deduced from the calibrations of spectral type- 
effective temperature or i n trinsic B-V color-effective tempera t ure given by 



Drilling fc Landolt 



mm 



Flower 



tool), 



de Jager fc Nieuwenhuijzen I ( 1l987l ) 



19961 ) and 



Popper 



Alonso et al. 



we estimated a temperature of 6 670±170 K for the 
frared colors - effective temperature calibrations of 



19801 ) is about 6 750±150 K. However 
jrimary sta r frorn the in- 



Tokunaga 



mm . Tem- 



perature uncertainty of the primary component results from considerations of 
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Figure 4: Sample of CCFs between T-LyrO-08070 and the RV template spectrum around 
the first and second quadrature. 

spectral type uncertainties, and calibration differences. As a weighted mean 
we adopted 6 700±150 K for the effective temperature of the primary star. 



3. Period determination 

We obtained five times of minimum light for the system during photo- 
metric observations. These times of minima are given in Ta ble 2 as averaged 
for the filters used. The times of minimum light obtained by [Fernandez et al. 
(120091 ) were also included in the list. A linear least-squares fit to the data 
yields the following ephemeris 

Mini = HJD 2455443.3232(2) + l'^. 18483992(41) x E 

where the bracketed quantity is the uncertainty in the last digit (s) of the 
preceding number. All uncertainties quoted in this work are standard errors. 
The residuals of the fit are plotted in Fig. 3 and show no indication of any 
form of period change in about five years. 



Table 2: Times of minimum light for T-LyrO-08070 and the residuals between the observed 
and calculated ti mes with the new ephem eris. 
References: (1) [Fernandez et al.l (|200g ): (2) This work. 



Time of minimum Cycle number O — C Reference 
(HJD 2500 000+) 



3996.6326±0.0001 


-1221.0 


-0.0010 


1 


3996.6333±0.0001 


-1221.0 


-0.0004 


1 


3996.6340±0.0001 


-1221.0 


0.0003 


1 


3996.6346±0.0001 


-1221.0 


0.0009 


1 


5443.3231±0.0004 


0.0 


0.0003 


2 


5444.5083±0.0002 


1.0 


0.0002 


2 


5455.1708±0.0004 


10.0 


-0.0008 


2 


5456.3550±0.0003 


11.0 


-0.0014 


2 


5716.4286±0.0003 


230.5 


-0.0002 


2 
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4. Analysis 

4-1. Radial velocity curve 

To derive the radial velocities for the components of binary system, the 
12 TFOSC spectra of the eclipsing binary were cross-correlated against the 
spectrum of G J 740, a single-lined MOV star, on an order-by-order basis 
using the FXCOR package in IRAF. The majority of the spectra showed two 
distinct cross-correlation peaks in the quadrature, one for each component 
of the binary Thus, both peaks were fit independently in the quadrature 
with a Gaussian profile to measure the velocity and errors of the individual 
components. If the two peaks appear blended, a double Gaussian was applied 
to the combined profile using de-blend function in the task. For each of the 
12 observations we determined a weighted-average radial velocity for each 
star from all orders without significant contamination by telluric absorption 
features. Here we used as weights the inverse of the variance of the radial 
velocity measurements in each order, as reported by FXCOR. 

We adopted a two — Gaussian fit algorithm to resolve cross-correlation 
peaks near the first and second quadratures when spectral lines are visible 
separately. Fig. 4 shows examples of cros-correlations obtained by using the 
largest FWHM at nearly first and second quadratures. The stronger peaks 
in each CCF correspond to the more luminous component which has a larger 
weight into the observed spectrum. 

The hehocentric RVs for the primary (Vp) and the secondary (Vg) com- 
ponents are listed in Table 3 and plotted in Fig. 5 against the orbital phase. 
The RVs listed in Tabic 3 arc the weighted averages of the values obtained 
from the cross- correlation of orders #4, #5, #6 and #7 of the target spec- 
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Table 3: Heliocentric radial velocities of T-LyrO-08070. 
The columns give the heliocentric Julian date, the or- 
bital phase (according to the ephemeris in Eq. 1), the 
radial velocities of the two components with the corre- 
sponding standard deviations. 



HJD 2400000+ 


Phase 


Star 1 


Star 2 








(7 


Vs 


(7 


55390.3850 


0.3204 


-70.5 


-3.5 






55391.4725 


0.2382 


-72.2 


-3.6 






55392.3905 


0.0130 


-34.3 


11.1 






55392.5751 


0.1688 


-73.0 


-3.7 






55393.2954 


0.7768 


13.8 


0.7 


-210.7 


9.9 


55393.3438 


0.8176 


9.6 


0.5 


-199.4 


11.1 


55394.5068 


0.7992 


12.9 


0.6 


-207.9 


12.8 


55396.3351 


0.3422 


-66.3 


-3.3 






55396.4813 


0.4656 


-40.2 


7.4 






55397.4438 


0.2780 


-73.9 


-3.7 


151.5 


10.4 


55398.5238 


0.1894 


-74.7 


6.2 






55398.5700 


0.2285 


-76.4 


-3.8 


149.9 


10.4 
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Figure 5: Radial velocity curve folded on a period of 1.1848 days. Points with error 
bars (error bars are masked by the symbol size in some cases) show the radial velocity 
measurements for the components of the system (primary: filled circles, secondary: open 
squares) . 



tra with the corresponding order of the standard star spectrum. The weight 
Wi = 1/crf has been given to each measurement. The standard errors of the 
weighted means have been calculated on the basis of the er rors (aj) in the R V 



values for each order according to the usual formula (e.g 



Topping 



(119721)) 



The (Tj values a re computed by FXCO R according to the fitted peak height, 



as described by 



Tonrv fc DavisI f ll979h 



First we analysed the radial velocities for the initial orbital parameters. 
We used the orbital period held fixed and computed the eccentricity of the 
orbit, systemic velocity and semi- amplitudes of the RVs. The results of 
the analysis follow: 7= -30±2 km s-\ Ki=45±3 and Js:2=190±12 
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Figure 6: The observed and computed light curves of T-LyrO-08070. 

km s~^with a circular orbit. Using these values we estimate the projected 
orbital semi-major axis and mass ratio as: asinz=5.50±0.06 Rq and q = 
^=0.237±0.008. 

4-2. Light curve modeling 

We used the most recent version of the ec lipsing binary light curve model- 



ing algorithm of 



Wilson fc Devinney 



( 1l97ll ) (with updates), as implemented 



in the PHOEBE code of Prsa & Zwitter (2005). The code needs some input 
parameters, which depend upon the physical properties of the component 
stars. In the light curve solution we fixed some parameters whose values can 
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be estimated from global stellar properties, such as effective temperature of 
the hotter component and mass-ratio of the system. The effective temper- 
ature of the primary star has already been derived from photometry and 
spectroscopy as 6 700 K and the mass-ratio as 0.237. A preliminary estimate 
for the effective temperature of the cooler component is made comparing the 
depths of the eclipses and mass-ratio. We adopted the linear limb-darkening 
coefficients from Van Hamme (1993) f or the primar y and secondary compo- 



nents, and the bolometric albedos from 



Lucy 



( 1l967l ). for each passband. The 



gravity brightening coefficients are taken as 0.32 for both components. The 
rotational velocities of the components are assumed to be synchronous with 
the orbital one. The adjustable parameters in the light curves fitting were 
the orbital inclination, the surface potentials, effective temperature of the 
secondary and the luminosity of the primary. Our final results of simultane- 
ous solution are listed in Table 4 and the computed light curves (continuous 
line) are compared with the observations in Fig. 6. The fractional radii (ri 2) 
given in Table 4 are the mean values. The proximity effects at out-off-eclipses 
are clearly seen. The uncertainities assigned to the adjusted parameters are 
the internal errors provided directly by the Wilson- Devinney code. 



5. Absolute dimensions and distance to the system 

The semi-amplitude of the radial velocities are obtained from the radial 
velocity analysis. The inclination of the orbit and the fractional radii of the 
components are determined from the light curve solutions. Therefore, we 
can compute the absolute parameters of the component stars. We present 
the fundamental physical parameters for the system T-LyrO-08070 with their 
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Table 4: Results of the simultaneous B, V, R, and I passbands light curve analysis for 
T-LyrO-08070. 



Parameters 


Adopted 


i° 




78.54±0.51 


Te/A (K) 


6 700 [Fix] 


Te/A (K) 


3 200±180 






3.650±0.012 


n2 




3.314±0.089 


ri 




0.285±0.004 


ri 




0.154±0.009 


Li 
(L1+L2) 

Li 
(L1+L2) 

Li 
(L1+L2) 

Li 
(L1+L2) 


(B) 
(V) 
(R) 

(I) 


0.998±0.008 
0.988±0.008 
0.910±0.008 
0.888±0.008 
0.006 



uncertainties in Table 5. The absolute parameters we obtained in th is study 



are significantly different from those obtained by [Fernandez et al.l ( 20091 ). 
because they derived the parameters for only the transit light curve and 
the rotational velocity of the primary star. In particular, the mass and 
ra dius of the secondary s tar are very different when compared with those 



of 



Fernandez et al. 



(120091 ) ■ They estimated a mass for the less massive star 
as 0.22, 0.25, 0.29 and 0.32 Mq , and radius as 0.25, 0.27, 0.29 and 0.31 
i?0from the stellar isochrones for the values of Fe/H=-1.0, -0.5, 0.0 and 
+0.5, respectively. Using orbit-rotation synchronization they find mass and 
radius for the secondary star as 0.240 M0and 0.265 Rq. Our radial velocity 
and light curve analyses yield a mass of 0.32 and a radius of 0.86 in solar 
units. 

In Fig. 7 we compare locations of the low-mass stars in the mass-radius 
diagram. The radius of the M3 dwarf is at least two times larger than 
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that predicted from the stellar models. The less massive components of 
the close binary systems having a mass smaller than the sun are generally 
oversized with respe ct to the theoretical models as indicated recently by 
Cakirli et al. I (120 lOf ). However, their effective temperatures are generally 
lower than those of single main sequence stars with similar mass. Larger radii 
and lower effective temperatures are attributed to the solar-like activities 
(see Morales, Ribas & Jordi (2008)). Large spots on the stars prevent energy 
output, therefore, the star expands. While the radius increases the surface 
temperature decreases in active stars in close binary systems. 

The BVJHK magnitudes of the binary and the comparison stars given in 
Table 1 allow us estimation of the color excess. Assuming the spectral type 
of the primary star as F3 we estimated the in terstellar reddening a s E(B - 
V)=0.22±0.02 mag using the tables given by briUing fc Landolt I (]200ol ). 
Using the radius and effective temperature of the primary star we compute 
its luminosity as 4.56 Lq. Then taking this value and bolometric absolute 
magnitude for the sun as 4.74 one may compute its absolute bolometric 
magi iitude as 3.08 "mag". The bolometric correction of -0.12 "mag" taken 



from 



Drilling fc Landolt I (|2000[ ) we find an absolute visual magnitude of 3.21 



"mag" for the more massive star. The apparent visual magnitude of 12.162 
extracting the light contribution of the secondary star to the total light, the 
absolute visual magnitude for the primary star and interstellar extinction 
of Ay =0.68 "mag" yield a distance to the system as 479±36 pc. However, 
JHK magnitudes lead to a distance of about 630 pc with larger standard 
deviation. This difference may be originated from the bolometric corrections 
taken from different tables (Girardi et al. (2002) and Kervella et al. (2004)). 
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Table 5: Fundamental parameters of T-LyrO-08070. 



Parameter 



Spectral Type 
a (Ro) 
(km s"^) 

in 

Mass (Mq) 
Radius (R©) 
log g (cgs) 

Teff (K) 



Primary Secondary 



F3±2 M3±l 
5.61±0.06 

-30±3 
78.5±0.5 
0.237±0.008 
1.37±0.23 0.32±0.04 
1.60±0.09 0.86±0.06 
4.17±0.03 4.08±0.05 
6700±150 3170±150 
{vsin i)obs (km s"^) 58±1 25±9 
{vsin i)caic. (km s'^) 69±2 27±2 
log (L/Lq) 0.67±0.06 -1.17±0.10 

d (pc) 479±36 
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Figure 7: Locations of the low-mass stars in the mass-radius diagram. The less massive 
companion of the T-LyrO-08070 is shown by a large square and the primary by dot with 
standard deviations. 
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